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Abstract 

We develop a method to calculate the prefactor in the expression for the 
bubble nucleation rate. A fermion with Yukawa coupling is considered where 
a step potential can be used as a good approximation in the thin wall limit. 
Corrections due to thicker walls are investigated by perturbing about the 
thin wall case. We derive the thermal one loop effective action, calculating it 
numerically, and find that the prefactor in the nucleation rate can both sup- 
press and enhance, for a given temperature, when the usual renormalisation 
conditions are applied to the effective potential. 
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I. INTRODUCTION 



Bubble nucleation can occur in a first order phase transition from tlie false to the true 
vacuum. The bubble nucleation rate per unit volume per unit time, written in the language 
of Coleman is r/\^ = Ae~^, where B is the classical Euclidean action of the bubble and 
A is the one loop contribution including zero and negative modes. 

The purpose of this paper is to present a method that enables one to calculate the 
nucleation rate in the thin wall limit. Techniques such as the derivative expansion break 
down for this type of background 0. The thin wall limit is a possible scenario in electroweak 
theory with multiple Higgs fields [0. It is also useful as a possible explanation for the 
generation of the baryon asymmetry we observe today 0. In the interests of brevity we 
focus on just the fermion fields. However, the general discussion includes scalar and spinor 
fields for completeness. The calculation of other fields, using the method presented, is 
currently under way 

Previous work on nucleation rates includes an analysis of the prefactor for scalar fields 
coupled to fermions by Gleiser, Marques and Ramos |^, who computed the determinant 
using the derivative expansion. Issues regarding which loop corrections should be included 
in the bounce solution were investigated. Garriga calculated the determinantal prefactor 
analytically, assuming the free energy can be approximated by that of a massless field living 
on the surface of a membrane. This was for a scalar field in the thin wall approximation 
(corrections due to thicker walls were also studied) at finite temperature, where enhancement 
of the nucleation rate was found. Kripfganz, Laser and Schmidt computed the prefactor 
for the one loop Higgs fluctuations at the electroweak phase transition, using the derivative 
expansion as an approximation. 

A particularly useful numerical method uses a theorem on functional determinants that 
can be found in Coleman's work [^]. Baacke and Kiselev develop an 'exact' numerical 
scheme to work out the one loop corrections for a scalar field at finite temperature, using 
Coleman's theorem. Thick and thin walls were considered but an infinitely thin wall {a = 
1 in their notation) was not. Hence, there is no direct method of comparison with 0. 
However, they also found an enhancement of the nucleation rate (depending on the choice of 
renormalisation scheme). Baacke [|10| then studied vector bosons in the 't Hooft-Feynman 
gauge, generalising Coleman's theorem to the matrix case. The results were compared to 
those in [0, where there was good agreement for small (thick wall) bubbles and strong 
deviations between the results for large (thin wall) bubbles, as expected. The nucleation 
rate was suppressed for this case. Also, Baacke and Siirig |1T| calculated the fermionic 
fluctuation determinant using Coleman's theorem and a gradient expansion for comparison. 
For an optimal choice of renormalisation, they found that the rate was enhanced at the 
electroweak temperature. 

Brahm and Lee []T2[ used a different numerical procedure based on the phase shift method 



and similar to the one we shall employ (using some old results due to Schwinger 1|T3[)- They 



computed the prefactor for scalar fields at finite temperature for the thin wall limit, based 
on the assumption that the surface free energy is equal to that of a domain wall (thick 
walls were also considered). The WKB approximation was used for the high energy modes 
to improve convergence of the exact result. The results were compared with the effective 
potential and derivative expansion approximations (see references therein). 
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Further work has been done by Miinster and Rotsch |T^, who calculated the prefactor 
in the nucleation rate for a scalar field using a Poschl- Teller potential and heat kernel coef- 
ficients to remove the ultraviolet divergences of the theory. Miinster, Strumia and Tetradis 



15| have compared this work to an entirely different method (not relying on saddle point 



evaluation) that uses coarse graining and renormalisation group techniques. They find good 
agreement between the two in the region of validity. Further work has been done using this 
coarse graining method and we refer the reader to [|T3| and the references therein. 

In our method we shall use phase shifts and relate the prefactor to the heat kernel. The 
theory is regulated by subtracting off the relevant heat kernel coefficients. We consider 
fermions with a Yukawa coupling and a step function profile to begin with. Using the step 
function gives an exact analytic expression for the phase shift, making the calculation more 
manageable. For the fermion case, the thin wall limit introduces problems of its own because 
the derivative of the mass term leads to a delta function in the potential. As far as we are 
aware, this is the first time this specific case has been investigated, although the method 
used in [jlO[ or |jl2[ could be applied. Results are at finite temperature (not only in the high 
temperature limit) and the method can be extended to the zero temperature case. Also, the 
technique involves a simple regularisation step, unlike methods based on Coleman's theorem 
which require evaluating uniform asymptotic expansions of the relevant field equations. 

The eigenvalues in the determinant are found using partial wave analysis and phase 
shifts, with the eigenmodes discretised by putting them in a sphere of large radius Q that 
we let tend to infinity (not to be confused with the bubble wall radius at R). Thus, we must 
impose boundary conditions on the fields. In the case of fermions, the correct eigenvalue 
problem requires mixed boundary conditions (see appendix A). 

The paper is organised as follows. In Sec. II A we relate the phase shift to the heat 
kernel and zeta function. In Sec. II B we calculate the thermal effective action leading 
to an expression for the prefactor in terms of the phase shift. In Sec. Ill we discuss how 
to consider corrections from thicker walls. In Sec. IV results are presented and in Sec. 
V we draw conclusions. In Sec. VI an appendix in three sections is given, with details 
on the calculation of the fermion phase shift, numerical zeta function regularisation and 
renormalisation respectively. 



II. HEAT KERNELS AND PHASE SHIFTS 

A. Nucleation &: regularisation 

We begin with the nucleation rates for the decay of a false vacuum at finite temperature 
due to an instanton (j) bubble- Any field that acquires a mass on the instanton background can 
contribute to the prefactor. In three dimensions [|l^] (see also comments in |^ & ||12||), 

det'[-V^ + m\<Pbume)] , . 

det[-W^ + m\(t)sym)] ■ 

Three zero eigenvalues, arising from breaking the Poincare symmetry, each contribute 

are omitted from the scalar determinant, as indicated 

by the prime. 



^=^1 n 
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In the thin wall limit we assume that (j) bubble has a bubble wall at some radius R, such 
that takes the false vacuum value at radii r > R and the true vacuum value at radii 
r < R, with a narrow transition region near r = R. For scalar bosons and fermions (and if 
calculated, the vector bosons), the relevant mass terms vanish in the false vacuum and are 
non-zero in the true vacuum, leading to a step function profile, 

f . ( m r < R /„x 
m r = < . (2) 

In the same limit, the scalar Higgs field masses differ little for large and small radii. This 
suggests that the Higgs contribution to the prefactor is smaller than the fermion contribution 
(and also any other fields). We will consider the accuracy of the thin- wall approximation 
later. 

The eigenvalues in the determinant can be found by using a partial wave analysis and 
phase shifts |T^. We first discretise the eigenmodes by putting them in a sphere of large 
radius Q. After separating the eigenmodes into radial functions and spherical harmonics, 
the radial parts asymptotically approach trigonometric functions of kr + 0, where is a 
constant phase depending on k and the angular momentum /. On the boundary, 

kn^ ^ mr — (f). (3) 

In the false vacuum, the potential is zero and we label the free eigenvalues k^^^ 

fcf fi^n7r-0W. (4) 

On letting Q —>■ oo (continuum limit), the above equations (^ and (|^) imply a relationship 
between the phase shift Si{k) = — 0^°-', the density of states gi{k) and gf\k), ||T3[| , 

for the radial modes. 

The difference between the heat kernels for the instanton and the true vacuum will be 

AK(t) = (e"'"* - e-'" (6) 

Using the density of states factor giik) and the degeneracy factor Xh 

AK{t) = r dke-^''Y.Xi{gi{k)-9f\k)). (7) 
Jo ; 

Substituting (|]) into the above equation and integrating by parts we obtain 

AK{t) = - dke-^'ktY^xMk), (8) 

where the degeneracy factor = (2/ + 1) in three dimensions. 

The heat kernel can now be used to regularise the determinants appearing in the prefactor 
A in the nucleation rate. We define the generalised C, function by 
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1 poo 

as) = —- f-HrKit)dt, (9) 
i [s) Jo 

The analytic continuation of ({s) then gives 

\ogA= J2 (±)AH^ (10) 

fields 

where we take ± for scalar and spinor fields respectively, and 

Aiy = iC'(o) + K(o)iog/^' (11) 

where fi is the renormalisation scale. 

For numerical work, the analytic continuation can best be performed by subtracting 
terms from the heat kernel. As t — 0, the heat kernel in d+1 dimensions has the asymptotic 
expansion [|I^ 

K{t) ~ t-('^+i)/2 J2 BX . (12) 

n=0 

The leading terms, which cause the poles in the C, function, can be removed by replacing the 
sum over phase shifts in equation (|^) by 

2^ XiOi - 2_, XiOi - ^^d+i^ T{^) ' ^ ' 

where the Bq coefficient cancels because it is equal to the free heat kernel K'^'^\ An infrared 
cutoff MjR must also be included, noting (see Appendix B) that the dependence on Mjji 
is illusory given that changing Mjji does not affect the value of /SW . The i?3/2 coefficient 
only occurs for fermions because squaring the Dirac equation leads to a delta function in the 
potential (for a step profile). Heat kernel coefficients have been calculated for distributional 
backgrounds pOf and we only quote the result below. 



For the step potential, standard expressions for the heat kernel coefficients give (& 



^1 = -:^J7T-^^^^. ' (14) 



^3/2 = ' (15) 



B2 = , ■ (16) 



2'^+i(rf + i)r 



For example, the phase shift [T^] for a scalar boson field is 

. ^ k Ji-i/2ikR) Ji+i/2{k'R) - k' Ji.i/2{k'R) Ji+y2{kR) 
'''' ' kNi_y,ikR)J,^,/2ik'R)-k'Ji.^/2ik'R)Ni^,/,ikR)' ^ ^ 



where k' = y/k"^ — w? and I = 0, 1, ... in three dimensions. The phase shift for fermions is a 
rather lengthy calculation that is left until Appendix A. 
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B. Thermal effective action 



Using the techniques of the last section we are now ready to calculate what is essentially 
the difference in the effective action for the true and false vacua AW. We refer the reader 
to for a detailed discussion of heat kernel methods at finite temperature for scalar and 
spinor fields. The thermal heat kernel K'^ can be expressed as an infinite sum of zero 
temperature heat kernels 

oo 

Rf^it I r, x; r', x') = ^ (±)" K{t \ r, x; r' + (3n, x') (18) 

n=— oo 

(where ± is for scalar and spinor fields respectively) and for an ultrastatic spacetime, the 
heat kernel can be factorised into temporal and spacial parts giving 

Kit I r, X] t\ x') = e-^^^K'-^\t \ x, x'). (19) 

v47rt 

It is then possible to show using the above relations that, 



oo 



i^^(t) = -^ ^ e-^K(^)(t). (20) 



n=— oo 



AW is related to the thermal heat kernel AK^{t) for scalar fields by 

1 /"OO rlf 

AW = -- / -trAK^'it). (21) 

2 Jo t 

Substituting (^) into the above equation, we obtain for scalars, 

AW = -^ ^ V e-^"'^- / dke-'"'kty2i2l + l)Siik). (22) 
2 ^0 t „^_^ V47rt TT Jo 

Then, using the fact that 

rt-h-^-'-'''dt = ^e--^^ (23) 
^0 k 



we get 



-n/3k 



(24) 



AW = -^ dkY^{2l + l)5i{k)-^ dkJ2i2l+l)5iik)J2 

2^^0 1=0 /=0 n=l 

The sum over n is standard, leading to the result 

l3 roo °° _ R roc °° iS,(k\ 

'^^y—li *i:(2;+i)^. (25) 

The first term is the zero point energy, that contains the ultraviolet divergences of the 
theory and hence 5i{k) (see equation ( pOD & (|T3D ). The second term is the temperature 
dependent part. The above expression can be derived using the density of states method 
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T^ . Thus, upon regulating the nonthermal part, using zeta function regularization and 



introducing a mass that we let tend to zero at the end of the calculation (see Appendix B) 
we have 



27r 



r dk ( ±{21 + m{k) - 2V^B,k - + logM'iR. (26) 



AW^^-^-fdkU2l,,)^^ 
n Jo eP'^ — 1 



(27) 



AW^(i/2)= / -«fv2)W- (28) 



where AW = AW^ + AW^ is the thermal effective action and we are working explicitly in 
three dimensions. (Note that AW is independent of Mfj^ and for the scalar boson, is 
zero.) 

For fermions the spinor effective action AW(^i/2) is related to the heat kernel AK'^^^2){'t) 
by (twice the scalar result for massive fermions and also a colour factor of three if the top 
quark is considered) 

f-oo dt 
Thus, 

AW^(i/2) = 4/5 / - E (-ire"-"ST^- / dke-'''ktY.m + l)Sf{k), (29) 

where 5/ = 5+ + 5_ (see appendix) and j = 1/2,3/2, ... in three dimensions (the factor of 
four comes from the trace over spinor indices). Applying (p3D gives 

Aft /-oo 00 S/9 /-oo 00 00 

Al^(i/2) = ^ / rfA; E 2(2j + 1) 6f{k) - ^ / rffc ^ 2(2j + 1) 5/(A:) ^(-1) 

^ ^0 ,=1/2 ^ i=l/2 n=l 

(30) 

Then, summing over n we have 

AW^(i/2) = f / dkY. m + 1) - f / rfA: E 2(2j + 1) (31) 

/I JU 3=1/2 ^ i=l/2 ~'~ 

Of course, we could have guessed this result from looking at (|^), taking into account 
the properties of spinors. It is fairly simple to derive the above equation using the density 
of states in the same way as in |0, using the spinor phase shift. Thus, 



= ^ dki E 2(2j + 1) 5f{k) - 2^B,k - 7r53/2 - 



logMfj,, (32) 



TT JO 



' — C^kt 2(2, + 1) -M^L. (33) 



i=l/2 " + ^ 



Note the opposite sign in the zero point energy (nonthermal) contribution to the one loop 
effective action as compared to the scalar case equation (^6]). 
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III. THICKER WALLS 



The phase shift method works equally well for any bubble profile, although a differential 
equation must be solved numerically to find the phase shift. In the general case, one could 
consider the difference between the phase shift and thin wall phase shift and then add this 
correction onto the effective action numerically. Alternatively, it is possible to estimate the 
corrections due to the non-zero thickness of the wall by perturbing about the thin wall case. 
For example, consider the scalar boson, for which we must solve 

- r-2(rV)' + m\ + ^-^^^u - k\ = -Vu, (34) 

where m is given by equation (|^), V is the correction due to a thicker wall and we define R 
such that r'^V{r)dr = 0. The Green's function is 

r(r r'] = -l kA-^Ui{r)u2ir') r<r' . . 

^ ' ' \kA-^U2{r)ui{r') r > r' ^ ^ 



and 



, ( ji{k'r) r <R ( , ( Cji{kr') + Dni{kr') r<R . 

\Aji{kr)-Bni{kr) r>R''^^^''> \ni{kr) r > R' ^ ' 



where k' = \Jk'^ — and jii^z) = y^z^'^ '^^^'^ Ji+{d-i)/2{z) in ci + 1 dimensions. 

One then imposes m oc as r — 5> and u oc Aji — B'rii as r — ^ oo. Then, the solution is 



u = Ui — / G{r,r')V{r')u{r') r'"^ dr' . (37) 
Jo 

Therefore, as r — > oo 

POO 

u^ui + kA-^U2 / V{r')ul{r') r'^ dr' . (38) 
Jo 

From this it is possible to show that the correction to the phase shift is 

k 

5W = -^^^y^ V{r)ul{r)r'dr. (39) 

Then, assuming that the Bessel functions change little as r varies over the bubble wall, they 
can be Taylor expanded about R, giving 



ji{k'R)j[{k'R) f V{r)r^dr, (40) 
Jo 



A^ + B"^ 

where the continuity of ui at the bubble wall has been used. For the scalar case, 

B = -kR\kji.^{kR)ji{k'R) - k'ji{kR)ji_^{k'R)) (41) 

and 

A = -kR\kni_,{kR)ji{k'R) - k'ni{kR)ji_,{k'R)). (42) 
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All the terms outside the integral in equation (1^) are numerical factors. For reasons dis- 
cussed in the conclusion we only add the correction onto the thermal part of the effective 
action, giving 

AiyW = [^Vhydr (43) 
R-^ Jo 

where a is given by 

1 /-oo - (2/ + 1) 2zz' . , 

If one assumes a ^tanh^ like potential m?'{r) for the thicker wall, then V = m?{r) — rn? ~ 
—m?/2 at the bubble wall radius R. Therefore, in terms of the width w of the bubble wall 
one obtains the approximate result AVT*^^) —al3m?w /2. We have calculated a numerically 
and found that it is of order 1 for a range of values of 13 / R. 

IV. RESULTS 

The thermal one loop effective action was calculated numerically for fermions. The phase 
shift is substituted into equations (|32D and (^31), where it is convenient to change variables 
k ^ z = kR for the step potential. Then the nonthermal part has only one parameter 
7] = m'j R^, where is the mass of the fermion. The thermal part has parameters rj and 
PfUf (given that (3/R = Prrif/y/rj) as independent variables. 

We work with equations (|3^) and (p3D using a numerical package. For each value of z 
the function is summed over I (or j), with I increasing up to a given L until the value of the 
function at z converges. The thermal part of the integral converges due to the exponential 
damping terms. When considering the nonthermal part of the function, we must check that 
the integrand has the correct k dependence after making the subtraction of the divergent 
quantities from the sum over the phase shift. This is a good check verifying that the heat 
kernel coefficients are correct. 

We integrate up to Z chosen to obtain the required accuracy (all results are accurate to 
1%). For large values of rj (and small (3m f), larger values of Z and L are needed to give 
convergence. The arctan function (from the phase shift) has problems with branches for 
large values of r] (whenever 6 hits ±7r), requiring numerical glueing of the phase shift. 

The nonthermal part of the fermion effective action can be written j3mfF{ri). Numer- 
ically, Firf) fits well to a power law dependence on rj (see Fig. P, giving, in the original 
variables, 

Aiy(f/2) = -1.51/3mji?2 + 0.32/3m)i?^ (45) 

where we have set the renormalisation scale /i = Mm = m/ as explained in Appendix C. 
The nonthermal part is plotted in Fig. The full one loop effective action for fermions 
plotted against rj is in Fig. |^ for various values of the parameter (3m f. 
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V. CONCLUSION 



We have presented a simple method to compute the prefactor in the expression for the 
bubble nucleation rate, applying this to infinitely thin walls. Analytic corrections due to 
thicker walls were considered by perturbing about the thin wall case. It is easy to see that 
these corrections do not effect the Bi heat kernel coefficient. Using the arguments from 
Appendix C, where it was shown that we only need to consider the renormalisation of the 
thin wall bubble, one can ignore the correction from B2. Thus, only adding corrections to 
the thermal part of the effective action should be a good approximation. 

In the context of the electroweak phase transition, we would also like to examine the 
vector bosons, using the method presented. This requires vector spherical harmonics and 
the relevant boundary conditions to calculate the phase shift. Then a full treatment of all 
the particle species at the phase transition can be worked out in the thin wall approximation, 
including analytic expressions for corrections due to thicker walls. 

The fermion contribution generally enhances the rate, but for large (3m f (when ^ = nif) 
it does not and becomes suppressive. In fact, at (3m j = 5.0 (see figure |^), the sign of the 
effective action changes for various values of 77 = mjR (for a fixed bubble wall radius). 
Choosing fi = mt and m/ = mt (the mass of the top quark) as in [|r^, then the log term 



cancels for the fermion determinant (see equation (^5])). Baacke and Siirig ||TT| found a 
negative contribution (enhance), whereas we find it can also be positive (suppress) for a 
temperature = 35GeV with (3mt = 5.0 and ^ 175GeV) lower than the that at the 
electroweak phase transition. However, there is no direct method of comparison between 



the two methods because in [|ri[] the thin wall limit was not considered. 

The renormalisation scale can be set by imposing conditions on the effective potential, 
calculated for constant background fields [|12|. In Appendix C it was shown that these 
conditions on the effective potential (that include the classical and regularised terms) lead 
to the choice = Mj^ = ^f- 

In the case of bubble nucleation, the bubble wall radius is found by extremising an action 
which includes the effective potential. The three dimensional action for the bounce solution 
is 

4:71 

B = ATx(3R^a (3R^e, (46) 

3 

where a and e are corrected surface and volume energy densities respectively. The nonther- 
mal part of the effective action (equation (145|) ) has a similar form and can be absorbed into 
a redefinition of a and e. The value of R should be presumably chosen to coincide with 
the extrema of the new action. Then, the nucleation rate is determined as a function of 
temperature by ^W'^ . 

The method described can also be used to calculate the one loop effective action in 
general, but in this case the Bq heat kernel coefficient must be renormalised away in the 
vacuum energy. The numerical procedure is simple and efficient. 
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VI. APPENDIX 

A. Fermion phase shift 

Here we present the calculation of the spinor phase shift. One must use the Dirac 
equation separated into radial and angular coordinates and also be careful in the way one 
works out the eigenvalues of the problem. Starting with 

{i'j ■ 6 — m)-?/'+ = Xip-, (47) 



{i'j ■ 6 + m)ip_ = -Xtp+, (48) 

so that upon squaring the Dirac equation we get the Klein Gordon equation. The radial 
components are then [|17| 



{E T m)g± + /4 + (j + 



(49) 



{E ± m)/± - + (j - 1)^ = ^A/^, (50) 



^± = (i't\.]-^ ^r = [: ; ), (51) 




where 



for j = I + 1/2. The solutions are spherical Bessel functions f± = A±jj_^_i and g± = C±jj_i 
(this can be seen easily by substituting ( ^D[ ) into (0)) where the constants are as yet 
undetermined. On substituting the power series for the Bessel functions into (|49|) and 
(^) one finds the relations 

{E T m)C± + A±k T XC^ = 0, (52) 



{E ± m)A± + C±k ± XA^ = 0. (53) 

The correct eigenvalue problem requires mixed boundary conditions [0 (also for a self- 
adjoint action) leading to /+ + 5^+ = /_ — = as r — > oo. It is also possible to show 
that there is a symmetry among the solutions such that /+ ^ /_ and (7+ ^ —9-- For the 
calculation we set E = and impose the boundary conditions with the above symmetries, 
noting that m ^ as r ^ 00. Using this information, it is possible to show that there are 
two solutions 

tan 5 = -i, (54) 
A+ ^ ^ 

for A_^. = ±C+ and = ±D_|_, where B± and C±. are the corresponding constants for nj_^i, 
the irregular solutions, that are spherical Neumann functions. 
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Then, matching the wavefunction at the junction and performing some algebra, we get 
the phase shifts, 

_ {k ± m) JjjkR) Jj+i{k'R) - k' Jjjk'R) Jj+i{kR) 
^ ~ {k ± m) Nj{kR) Jj+i{k'R) - k' Jjik'R) Nj+^{kR) ' ^ 

where all symbols have their usual meanings as in the rest of the paper. 

The j = 1 — 1/2 modes have symmetry such that f± g± changes the boundary 
conditions and equations into the j = 1 + 1/2 case. Thus we have the above phase shifts 
with a total degeneracy 2(2j + 1), where j = 1/2,3/2, .... All results can be generalised to 
the four dimensional case. 



B. Numerical zeta function regularisation 

Consider the finite temperature case. The nonthermal part of the heat kernel is 

AKoit) = -^=- r dkke-^'"+'''"^'T6Jk)t, (56) 

where all symbols have the same conventions as in the rest of the paper and we have intro- 
duced a mass term M. The zeta function is then 

ACo(.) = ^-^^^^ rdkk{e + Mr^~'^Y.^.ik). (57) 
^Ji4:^T)^ r(s) Jo ^ 

The integral converges for s > 1. In order to apply the zeta function to numerical work, we 
subtract terms from the asymptotic expansion of the heat kernel and define 

P{S) = - /^^^'^}^ r + ^')~^'^ IE ^7 - ^V^B.k - 7ri?3/2 - v^^l . 

TT ./(47r) r(s) Jo ' k } 



(58) 

Performing the integration over k, and adding back the subtracted terms, one finds 

AC„(.) ^ Pis) + ^M--j^ + JLlilZllM^-^-B,. + ^M-B. (69) 

This is a regular function at s = 0. Taking the derivative with respect to s (for s = 0) and 
the limit M 0, we obtain 

Ac;(o) ^'Ij-^ki^p^- 2V^B.k - - n0=lj - ^ m 



where Mm is an infrared cutoff. It is important to note that if one subtracts the difference 
between ACq(O) for different values of infrared cutoff {Ma and Mb say) and using the fact 
that 



oo 



dk I , ^ - , ^ I = log^, (61) 



it is clear that this difference is zero. Thus, equation (|60|) is independent of M/r. 
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C. Renormalisation 



The renormalisation scale /i represents an arbitrary parameter that can be expressed 
in terms of other constants by imposing conditions on the effective potential (see for 
example). We shall show that the usual renormalisation conditions for the effective potential, 

V'{(t)^,ac) = 0, 

^"(0.ac) = m,, (62) 

where (f)yac is the classical vacuum expectation value and nih is the mass of the Higgs field, 
give the optimal choice, removing the logarithmic term. 

The first two terms in the zero point energy (|45D come from the finite, subtracted, part 
of the effective action (see equation (|58|)). A similar approach was used in where the 



divergent effective action was regularised by subtracting terms proportional to 1st and 2nd 
powers of the background potential, leaving a finite contribution to the fluctuation deter- 
minant. The two divergent terms were then regulated by using a momentum cutoff and 
renormalised by applying conditions on the propagators. This procedure is equivalent to the 
subtraction of the heat kernel coefficients up to and including B2 (see equation ([58|)). How- 
ever, zeta function regularisation is unusual, as compared with most regularisation schemes, 
in that it does not need infinite counter terms. Thus, in equation (|60D only the last term, 
containing the logarithm, requires renormalisation. 
The effective potential V = Vq + Vi will be defined 



V = / /3d^x 



\{ml + bml)c\? - (^ + bg)c\? + ^(A + 5A)0^ 



'log— ^, (63) 



where the terms in the square brackets make up the classical potential (including countert- 
erms), the last term is the (^-regularised 1-loop contribution with 

= ACo(O) = J d'x ^ m)0^(f ), (64) 



and we have included the renormalisation scale /i (see equation (|Tl])). 



In our case the renormalisation must be applied to a step function field configuration, 
but the ultraviolet divergences are independent of the background profile. This is intuitively 
obvious because k ^ 00 equates to the short distance behaviour of the theory. Thus, the 
B2 heat kernel coefficient (cx 4>'^{x)) need only be considered for (j){x) = constant, i.e. the 
effective potential. 

By applying the renormalisation conditions (|62|) to the effective potential, it is simple to 
show that the counter terms are given by 6m1 = 6g = and 

a = ^log#. (65) 

On substituting the above counterterms into the effective potential it is clear that the 
logarithmic term will cancel. Naturally, for our numerical procedure we can choose Mjpt = 
nif because it was shown in Appendix B that the infrared cutoff was independent of the 
result (equation (|60|)). 



13 



With our choice of renormahsation we can compare with a similar work [0, where the 
scalar Higgs field was investigated. In the case of the thin wall limit they only considered the 
surface free energy. However, comparing with our zero point surface energy contribution, 
for a scalar boson [^, we find it is of a similar order of magnitude and also suppresses the 



nucleation rate. This contribution is considerably smaller than that made by a fermion with 
a Yukawa coupling which, generally, enhances the nucleation rate. 
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FIG. 2. The thermal part of the one loop contribution to the fermion effective action AW^|^ 
plotted against 77 for values of /?m/ (from bottom to top); 0.5, 1.0 & 5.0. 
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FIG. 3. The full one loop contribution to the fermion effective action plotted against 

r) for values of /3m/ (from bottom to top); 0.5, 1.0 &; 5.0. 
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